AdOO  31  IZT'llSH 


ESD-TR-66-111 


MTR-92 


ESD  RECORD  COPY 

RETURN  TO 

SCIENTIFIC  &  TECHNICAL  INFORMATION  DIVISION 
£SIi).  BUILDING  1211 


ESD  7 

ESTI  Cali  No 
Copy  No. 


At  5. 


«-4-3PST 


cys. 


SPHERICAL  ZONE  BOUNDARIES  IN 
PERSPECTIVE  PROJECTION 


OCTOBER  1966 


B.A.  Forest 

Prepared  for 

DEPUTY  FOR  ENGINEERING  AND  TECHNOLOGY 

DIRECTORATE  OF  COMPUTERS 
ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
L.  G.  Hanscom  Field,  Bedford,  Massachusetts 


Distribution  of  this  document  is  unlimited. 


Project  7070 

Prepared  by 

THE  MITRE  CORPORATION 
Bedford,  Massachusetts 

Contract  AF  19(628)-5 165 

i\hb  M 


This  document  may  be  reproduced  to  satisfy  official 
needs  of  U.S.  Government  agencies.  No  other  repro¬ 
duction  authorized  except  with  permission  of  Hq. 
Electronic  Systems  Division,  ATTN:  ESTI. 


When  US  Government  drawings,  specifications,  or 
other  data  are  used  for  any  purpose  other  than  a 
definitely  related  government  procurement  operation, 
the  government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever;  and  the  fact  that  the 
government  may  have  formulated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications, 
or  other  data  is  not  to  be  regarded  by  implication 
or  otherwise,  as  in  any  manner  licensing  the  holder 
or  any  other  person  or  corporation,  or  conveying 
any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  invention  that  may  in  any  way  be  related 
thereto. 


Do  not  return  this  copy.  Retain  or  destroy. 


ESD-TR-66-111 


MTR-92 


SPHERICAL  ZONE  BOUNDARIES  IN 
PERSPECTIVE  PROJECTION 


OCTOBER  1966 


B.A.  Forest 


Prepared  for 

DEPUTY  FOR  ENGINEERING  AND  TECHNOLOGY 

DIRECTORATE  OF  COMPUTERS 
ELECTRONIC  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
L.  G.  Hansoom  Field,  Bedford,  Massachusetts 


Distribution  o(  this  document  is  unlimited. 


Project  7070 

Prepared  by 

THE  MITRE  CORPORATION 
Bedford,  Massachusetts 

Contract  AF  19(62S)-5 1 05 


FOREWORD 


The  author  thanks  Dr.  S.  Okada  for  the  interest  he  has  shown  in  this 
report  and  the  many  useful  suggestions  he  has  made. 


REVIEW  AND  APPROVAL 


This  technical  report  has  been  reviewed  and  is  approved. 


rv^i 


GEORGE  E.  VRANESH 
1st  Lt.  USAF 


7070  Project  Officer 


ii 


ABSTRACT 


This  report  considers  the  projection  of  spherical  zone  boundaries  onto 
an  image  plane  using  the  "perspective  projection".  Special  cases  of  this  type 
of  mapping  include:  the  stereographic,  the  gnomonic,  and  the  orthographic 
projections. 

In  addition  to  proving  that  the  image  of  any  spherical  zone  boundary  is 
a  conic  section,  the  formulation  of  the  parameters  of  that  conic  section  is  given. 

An  important  aspect  of  the  projection  technique  lies  in  the  capability  of 
having  any  point  of  the  sphere  as  the  center  of  projection  and  any  great  circle 
through  that  point  as  the  centerline  of  the  mapping. 
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SECTION  I 


INTRODUCTION 

In  the  development  of  analytical  or  visual  aids  for  aerospace  studies 
there  has  been  a  growing  need  for  a  technique  to  produce  maps,  geographic  re¬ 
ference  systems,  orbital  track  overlays,  etc. ,  quickly  and  accurately. 

This  report  presents  a  discussion  of  mapping  called  "perspective  projec¬ 
tion",  and  its  application  to  the  projection  of  spherical  zone  boundaries  onto  a 
plane. 

In  Section  II  basic  projection  formulas  are  derived,  and  the  special  cases 
related  to  well  known  map  projections  are  given. 

In  Section  HI  the  application  to  the  earth  is  discussed,  and  Section  IV 
gives  the  technique  for  making  any  given  geographic  location  the  center  of  pro¬ 
jection. 

Section  V  develops  the  general  theory  of  projecting  spherical  zone  bound¬ 
aries.  In  addition  to  proving  that  the  image  is  always  a  conic  section  (allowing 
for  degenerate  cases),  the  conditions  under  which  certain  conic  sections  arise 
are  given  in  Section  VI. 

In  Section  VII  each  conic  section  is  discussed  in  detail.  The  formulation 
is  given  in  such  a  way  that  programming  logic  will  easily  follow. 

It  is  intended  that  the  perspective  projection  will  be  set  up  in  a  master 
program  to  produce  maps,  geographic  reference  grids,  station  reference  grids, 
and  various  curves  for  use  as  overlays,  such  as  orbital  tracks,  visibility  zones, 
satellite  eclipse  zones,  etc.  It  is  intended  that  a  Benson-Lehner  plotter  will  be 
used. 
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Many  of  the  programs  mentioned  above  have  been  written  and  are  opera¬ 
ting.  These  will  be  discussed  in  future  working  reports. 
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SECTION  II 


THE  PROJECTION 


DESCRIPTION 

In  the  general  case  one  has  an  object,  O,  in  space;  a  point  of  projec¬ 
tion,  P;  and  an  image  plane,  I.  (See  Figure  1. )  A  line  is  drawn  to  a  point  T 
on  the  object  (or  preimage).  If  the  line  PT  (extended)  intersects  the  image 
plane  I,  the  point  of  intersection  T'  will  be  called  the  image  of  T. 


The  shape  of  the  image  depends  upon  the  shape  of  the  object,  the  orien¬ 
tation  of  the  image  plane  with  respect  to  the  object,  and  the  position  of  the  point 
of  projection. 
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LIMITATIONS 


This  discussion  is  restricted  in  the  following  ways: 

(a)  only  those  object  curves  in  space  are  considered  which  are 
generated  by  the  intersection  of  a  plane  with  a  given  sphere 
and 

(b)  the  image  plane  is  perpendicular  to  the  line  joining  the  point 
of  projection  and  the  center  of  the  given  sphere. 

FORMULATIONS 

The  geometric  setup  in  space  is  formulated  as  follows  (see  Figures 
2  and  3). 
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Figure  2.  Projection  Related  to  a 
Unit  Sphere 
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Figure  3.  Cross  Section  of 
Figure  2 


The  usual  (x,  y,  z)  coordinate  system  is  used  with  unit  vectors  1,  j,  1c 
giving  the  positive  direction  of  the  x,  y,  and  z  axes  respectively.  This  notation 
is  used  with  respect  to  the  object.  Points  of  the  image  will  be  designated,  in 
the  same  frame,  using  coordinates  x'  ,  y' ,  and  z' . 

A  unit  sphere  with  center  at  the  origin  is  used  whose  equation  is 

2  2  2 

x  +  y  +  z  =  1.  (1) 

The  point  of  projection  p  is  given  by 

P  =  P  •  £  •  (2) 

The  image  plane  is  given  by 

z  =  w.  (3) 

A  point  T  of  the  object  in  space  has  coordinates  (x,  y,  z),  and  a 
position  vector 

1'  =  xi  +  yj  +  z£  ,  (4) 

or,  in  spherical  coordinates  (a,  /3 ) 

T  =  cosa  cosi3i  +  sina  cos/3j  +  sin  /3k  (5) 

Figure  3  shows  the  geometry  resulting  from  passing  a  plane  through 
0,  the  center  of  the  sphere,  the  point  of  projection,  p,  and  the  point  T 
to  be  projected.  The  position  vector  of  T’  ,  the  image  of  T,  is  given  by 

T’  =  p£  +  [T  "  p£]  for  P  +  t3  (6) 

where  T^  ,  T^  and  T^  are  the  x,  y,  z  components  of 

Primed  letters  shall  be  used  to  indicate  points  in  the  image  plane  cor¬ 
responding  to  points  of  the  sphere. 
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In  the  image  plane  I,  we  find  the  coordinates  (x' ,  y ' )  of  T'  as, 


yf  =  t'  •  j  =  ‘  t2  p^t3  (?) 

3 

Thus  any  point  of  the  sphere  (for  which  p  4  T  )  may  be  projected  onto  the 

O 

image  plane  from  point  P  by  using  Equations  7. 

PARAMETERS  P  AND  W  RELATED  TO  STANDARD  MAP  PROJECTIONS 

Corresponding  to  certain  values  of  p  and  w  are  standard  map  pro¬ 
jections  as  defined  by  cartographers.  The  choices  of  p  and  w  and  the  re¬ 
lated  map  projections  are  listed  in  Table  I. 

Details  of  a  specific  projection  may  be  obtained  from  Reference  1. 

Table  I 

Parameters  P  and  W  Related  to  Standard  Map  Projections 


p 

w 

Projection 

Scale 

-1 

0 

Stereographic 

1 

unit  =  3963 

mi. 

0 

1 

Gnomonic 

1 

unit  =  3963 

mi. 

go 

0 

Orthographic 

1 

unit  =  3963 

mi. 

p 

w  =  l/p 

Geometric  projection 

1 

unit  =  3963 

mi. 

Figure  4  shows  the  geometry  of  the  projection  for  each  case,  using  a 
setup  as  indicated  in  FORMULATIONS,  beginning  on  page  4. 
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STEREOGRAPHIC  PROJECTION 


Z 


r 


pV  s"t' 

(  b) 


GNOMONIC  PROJECTION 


Z 


(d) 

GEOMETRIC  PROJECTION 


Figure  4.  Geometrical  Setup  for  Various  Values  of  P  and  W 
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SECTION  m 


APPLICATION  TO  THE  EARTH 

If  the  unit  sphere  as  given  by  Equation  (1)  represents  the  earth' s  surface, 
the  following  is  defined: 

(a)  the  north  pole  normally  has  position  vector  lc  and  the  south 
pole  -6. 

(b)  the  Greenwich  meridian  lies  in  the  x-z  plane;  the  upper  meridian 
of  Greenwich  is  that  one-half  of  the  great  circle  which  contains 
Greenwich. 

(c)  longitude  (east)  and  latitude  will  correspond  to  the  usual  spherical- 
coordinate  angle  measurements. 

A  point  T  of  the  sphere,  then,  has  a  position  vector, 

*  A  ~  ^ 

T  =  cosp  cosM  +  cosp  sinAj  +  sinpk 

where 

a  =  east  longitude  and  P  =  latitude  of  the  position. 

Using  Equations  (7)  it  is  clear  that  if  p  4  1,  the  north  pole  projects 
always  to  the  point  (0,  0)  in  the  image  plane  and  the  projection  is  referred 
to  as  the  "normal"  projection. 
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SECTION  IV 


CHANGING  THE  CENTER  OF  PROJECTION 

The  foregoing  formulation  restricts  one  to  a  map  in  which  the  North  Pole 
is  the  center  of  projection  and  the  Greenwich  meridian  is  the  centerline  of  the 
map.  It  is  desirable  to  have  other  geographic  locations  as  the  center  of  the  pro¬ 
jection  and  also  to  have  any  great  circle  through  this  point  as  the  centerline  of 
the  map.  This  section  sets  forth  the  theory  and  formulation  necessary  to  accom¬ 
plish  these  results. 

As  depicted  in  Figure  5,  the  new  center  of  projection  is  chosen  as  the 
point  C^A where  ^A^,  are  respectively  the  east  longitude  and  lati¬ 
tude  of  C.  A  certain  great  circle  through  C  is  determined  by  choosing  a 
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second  point  B^A^,  not  at  C.  The  great  circle  containing  arc  CB 
makes  an  angle,  a,  with  the  meridian  of  C.  Note  here  that  a  is  the  azimuthal 
angle  of  B  measured  clockwise  from  the  meridian  of  C.  Clearly  0  ^  a  ^  2n. 
Finally,  the  Greenwich  meridian  contains  arc  GP,  and  the  coordinate  axes  x 
and  z  contain  respectively,  the  segments  OG  and  OP  with  O  as  the  center  of 
the  sphere.  The  axis  of  y  makes  a  right-handed  system  with  the  x  and  z  axes. 

A  transformation  R  is  now  developed  which  will  do  the  following,  in 

order. 


It  will  rotate  the  sphere  about  the  z-axis  so  that  point  C  (or  the  meri¬ 
dian  of  C)  moves  through  the  angle,  -A  ;  C  then  lies  on  the  prime  meridian. 

L/ 

This  first  rotation  is  given  by  the  matrix, 


It  will  then  rotate  the  sphere  about  the  y-axis  through  the  angle, 

-  ^r/2  -  bringing  the  point  C  to  the  position,  (0,  0,  1),  formerly 

occupied  by  the  North  Pole.  This  second  rotation  is  given  by  the  matrix, 


(bos^r/2  -  <t>^j  0 

0  1 

sin(7r/2  -C)  0 


sin(7r/2  -  C)' 
0 

cos(7r/2  -  C) 


sin  bc 
0 

cos  0 c 


0 

1 

0 


-  cos 
0 

sin 


Finally,  there  is  a  rotation  of  the  sphere  about  the  z-axis  through  the 
angle  (7r  +  o)  carrying  the  point  B  counterclockwise  to  a  position  on  the 
prime  meridian.  This  last  rotation  is  given  by  the  matrix, 
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43os(7T  +  a)  -  sin(7T  +  a)  0 
sin(?r  +  a) 

0 


+  sin  a 
-  cos  ca 
0 


The  angle  a  is  found  as  follows. 

Let  6  and  B  be  the  position  vectors  of  points  C  and 
spectively  and  P  =  K  be  the  position  vector  of  the  pole,  then 


re- 


A  A  A  A 

(C  x  P)  .  (C  x  B) 

COS  OL  =  - - - - - - 

cos  P  .  sin  6 

u 

cos  0  sin  0  -  cos  0  sin  0  cos  A  -  A 

U  ±i  B  C  C  B 

sin  6 


__  \  A  A 

where  6  =  cos  (C  •  B)  is  the  distance  from  C  to  B  in  radians.  This 
formula  makes  use  of  the  fact  that  the  azimuth  angle  measures  the  dihedral  angle 
between  the  planes  containing  arcs  CB  and  CP.  This  dihedral  angle  is  the 
angle  between  vectors  drawn  perpendicular  to  the  planes.  (One  must  be  careful 
of  the  direction  of  these  vectors. ) 


From  the  spherical  triangle  CPB 

sin/A^  -  A 
sin  a  =  - 


we  get  the  relationship 

c)  eos  '’b 


sin  6 


(10) 


There  is  now  sufficient  information  to  assign  quadrants. 
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The  product  matrix  R 


BCD  follows  as: 


/cL  cl  3.  ^ 

11  12  13 


R  =1  a21  a22  a23 


a31  a32  a33> 


and  in  which 


all 

^cos 

a  sin 

00  cos 

X 

c 

+ 

sin 

a 

sin 

xc)' 

ai2 

=  - 

cos 

a  sin 

00  cos 

X 

c 

+ 

sin 

a 

cos 

Ac’ 

ai3 

= 

cos 

a  cos 

v 

a21 

=  __ 

sin 

a  sin 

00  COS 

X 

c 

+ 

cos 

a 

sin 

V 

a22 

=  - 

^sin 

a  sin 

00  sin 

xc 

+ 

cos 

a 

cos 

xc> 

a23 

= 

sin 

a  cos 

0 

c, 

a31 

= 

cos 

0q  cos 

x  , 
c’ 

a32 

= 

cos 

00  sin  A  ,  and 

sin  P, 


33 

Thus,  any  position  vector  of  a  given  point  A  of  the  sphere  will 
rotated  into  a  new  position  vector 


(11) 
be 


A  =  rAc. 


The  position  vector  A^X^,  0q)  before  rotation  is  given  by 
'  cos  X  f  +  cos  0  sin  X  f  +  sin  0„  6. 

0  0  oo  0 


A 

Aq  =  cos  P 
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Applying  the  transformation  given  by  Equations  (11)  the  components 
(a,  b,  c)  of  A  are  derived  as: 


a  = 

an 

COS 

% 

cos 

\  + 
0 

ai2 

cos 

0o 

sin  Ao 

+ 

ai3 

sin 

b  = 

a21 

cos 

0o 

cos 

Ao  + 

a22 

cos 

0O 

sin  A0 

+ 

"23 

sin 

and 

c  = 

a31 

cos 

0o 

cos 

A  + 

0 

a32 

cos 

0» 

sin  A0 

+ 

a33 

sin 

0o- 

Finally,  on  applying  Equations  7 


The  point  (x'  y' )  is  the  image  of  the  point  A. 


for  p  ^  c 
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SECTION  V 


PROJECTING  SPHERICAL  ZONE  BOUNDARIES 


In  this  section  the  projection  of  spherical  zone  boundaries  from  the 
sphere  onto  the  image  plane  is  considered.  A  spherical  zone  boundary  is  a 
circle  in  space  and  is  the  intersection  of  a  plane  with  the  sphere.  The  plane 
is  uniquely  determined  by  a  unit  vector 

A  =  ai  +  bj  +  clc 


normal  to  the  plane  and  a  fixed  point  B  in  the  plane  whose  position  vector  is 

B. 

If  1c  is  the  position  vector  of  any  point  (x,  y,  z)  in  the  plane,  then 

(x  -  3)  •  A  =  0 
defines  the  plane  (see  Figure  6). 


Figure  6. 


Geometry  of  Defining  the  Plane 
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Equation  (12)  reduces  to 

A  •  x  =  B  •  A  =  d 

where  d  is  the  directed  distance  from  the  origin  to  the  plane. 

2 

It  is  required  that  -  1  <  d  <  1.  If  d  =1,  then  the  plane  is 
tangent  to  the  sphere  and  we  have  a  null  circle.  If  d  =  0,  the  plane  contains 
the  center  of  the  sphere  and  the  zone  boundary  is  a  great  circle.  In  fact,  the 
zone  boundary  is  a  great  circle  if  and  only  if  d  =  0. 

Assume  that  the  point  to  be  projected  lies  on  the  plane  given  by  Equa¬ 
tion  (12)  viz: 

ax  +  by  +  cz  =  A‘B  =  d  (13) 

This  plane  will  be  referred  to  as  the  generating  plane,  G. 

Now,  Equations  (7)  relate  a  point  T(x,  y,  z)  in  space  to  a  point 
T'  (x'  y' )  in  the  image  plane.  If  the  point  T  lies  in  the  plane  given  by 
Equation  (13) , 

(PC-d)-x'  , 

[  c(p  -  w)  -  (ax'  +  by' )] 

_  _ (pc  -  d)  y'  ,  and 

^  [  c(p  -  w)  -  (ax'  +  by' )] 

„  _  d(p  -  w)  -p(ax'  +  by' )  . 

c(p  -  w)  -  (ax'  +  by' )  ' 
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If  the  point  T(x,  y,  z)  is  further  constrained  to  lie  on  the  sphere 
given  by  Equation  (1) ,  the  following  condition  exists, 


x,2[(d-pc)2  +  a2(p2-l)]  +  y'2[(pc-d)2  +  b2(p2-l)] 
+  2ab(p2-l)  x'  y'  +  2(p-w)  (c-pd)  (ax'  +  by' ) 


(p-w) 


•2(d2-c2)  = 


(15) 


Equation  (15)  clearly  defines  the  image  to  be  a  conic  section  of  the  form 


Au  +  Buv  +Cv  +  Du+Ev+F  =  0. 


(16) 


By  using  the  methods  set  forth  in  Reference  2,  page  256,  and  the  no¬ 
tation  of  Equation  (16),  we  find  the  following  equation: 

B2  -  4AC  =  -  4  [  (d  -  pc)2]  [(d  -  pc)2  +  (a2  +  b2)  (p2  -  l)  ]  (17) 


and  the  discriminant  (see  the  Appendix). 
2A  B  D 


A  =  — 


B  2C  E 


D 


2F 


=  -  4(p  -  w)  (d  -  pc) 


(i-A 


(18) 


16 


The  type  of  conic  that  one  gets  for  an  image  is  summarized  in  Table  II. 


Table  H 

Type  of  Conic  Sections  Related  to  Values  of  -  4AC  and  A 


2  .  .  _ 

B  -  4AC  <  0 

B2  -  4AC  =  0 

B2  -  4AC  >  0 

A  ^  0 

Ellipse  if  AA  <  0 

No  locus  if  AA  >  0 

Parabola 

Hyperbola 

A  =  0 

Point  (null  ellipse) 

2  parallel 
lines,  a 
single  line, 
or  no  locus 

2  intersect¬ 
ing  lines. 

Now  analyze  Equations  (17)  and  (18)  to  determine  the  geometric  situa- 

o 

tion  for  which  A  =  0  and  B  -  4AC  is  positive,  negative,  or  zero.  For  this 

and  later  applications,  it  is  convenient  to  write  Equation  (15)  in  a  different  form, 

viz. : 

(d  -  pc)2  (*'  2  +  y' 2)  +  [  (P  +  1)  (ax'  +  by' )  -  (d  +  c)  (p  -  w)] 

[  (p  -  1)  (ax'  +  by' )  -  (d  -  c)  (p  -  w)]  =  0.  (19) 


To  simplify  the  analysis  and  the  geometry,  it  will  be  helpful  to  rotate 
the  sphere  about  the  z-axis  in  such  a  way  that  the  vector  A  goes  into  a  vector 
A'  which  lies  in  the  x-z  plane  and  A'  •  i  >  0. 


Using  the  matrix 


R  = 


where  m  =  +  b^ 


(20) 
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we  have 


A'  =  RA  =  mi  +  ck 


(21) 


Consider  first,  the  case  A1  •  i  ^  0  ,  that  is  the  generating  plane  is 
not  parallel  to  the  x-y  plane. 

Figure  7  shows  the  trace  in  the  x-z  plane.  The  line  GR  is  the  trace 
of  the  generating  plane,  G  .  QQ'  is  the  projection  of  that  spherical  zone 
boundary  generated  by  the  intersection  of  plane  G  with  the  sphere.  I  is  the 
trace  of  the  image  plane,  and  p  is  the  point  of  projection. 

A  A  A  A  A 

We  shall  have  use  for  a  unit  vector  N  such  that  N  •  A  =  0;  N  •  j  =  0; 
the  desired  vector  is 


A  A  A 

N  =  -  ci  +  mk 


in  the  rotated  system,  and 


-ac  7  be  7 


N  =  -  i  —  j  +  mk 

m  m 


(22) 


in  the  nonrotated  system.  Clearly  the  vector  N'  is  normal  to  A'  ,  lies  in 
the  x-z  plane,  and  has  a  positive  z-component  (m  >  0). 

It  is  clear  that  the  points  Q  and  Q'  are  the  points  on  the  circular 
zone  boundary  which  have  the  maximum  and  minimum  z  coordinates  respec¬ 


tively.  (Note  that  since  A'  is  normal  to  plane  G,  the  rotation  brings  the 
plane,  G,  perpendicular  to  the  x-z  plane,  hence  the  trace  GR  is  also  the 
projection  of  plane  G  on  the  x-z  plane. ) 

The  position  vectors  of  Q  and  Q'  are, 


and 


0Q'  =  Q 


(23) 
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Since  d  ^  1,  A  is  never  positive.  Consider  now  each  of  the  factors 
in  the  right  member  of  Equation  (18): 

(a)  p  -  w  =  0  if  p  =  w; 

that  is,  if  the  image  plane  contains  the  point  of  projection.  In  a  practical 
situation  p  #  w; 

(b)  1  -  d2  =  0  if  d  =  ±  l; 

that  is,  if  the  generating  plane  is  tangent  to  the  sphere.  In  this  case  the 
degeneracy  arises  from  the  fact  that  a  point  is  projected.  In  particular,  the 

A 

point  is  either  that  whose  position  vector  is  A  for  d  =  +  1  or  the  antipode 
of  this  point  for  d  =  -  1  . 

(c)  d  -  pc  =  0 

if  the  generating  plane  contains  the  point  of  projection.  This  is  shown  as 
follows  (see  Figure  7).  The  equation  of  the  line  RG  is 

rax  +  cz  =  d  .  (24) 

The  z-intercept  of  this  line  is,  in  general,  the  point  at  which  the  generating 
plane  cuts  the  z-axis.  Thus  (see  Figure  7). 

— ►  d  - 

OR  =  —  k  for  c  +  0  . 
c 

If  the  generating  plane  contains  the  point  P,  we  have  d/c  =  p  or 

d  -  pc  =  0  .  (25) 

Thus  d  -  pc  =  0  if  and  only  if  the  generating  plane  contains  the  point  of 
projection. 

If  c  =  0,  the  generating  plane  is  parallel  to  the  z-axis  and  from 
Equation  (25)  we  infer  that  it  is  necessary  that  d  =  0  in  order  for  the 
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generating  plane  to  contain  the  point  of  projection.  In  summary,  a  degenerate 
case  is  always  present  when  a  great  circle  passes  through  the  center  of 
projection. 

Consider  now  the  factors  of  the  right  member  of  Equation  (17): 

2 

(a)  (d  -  pc)  2  o.  We  have  considered  this  before  as  to  the  equality 
with  zero.  Since  this  factor  is  never  negative,  it  will  not  affect  the  sign  of 
(B2  -  4AC). 

(b)  F  =  [(d  -  pc)2  +  (a“  +  b2)  (p2  -  l)] .  This  factor  is  not  so 

2 

easily  disposed  of .  If  (d  -  pc)  #  0,  then  B  -  4AC  is  positive,  negative, 
or  zero,  according  as  F  is  negative,  positive,  or  zero,  respectively. 

2 

If  p  si,  then  F  ^  0.  Thus  only  those  values  of  p  for  which  the 

2 

point  of  projection  lies  inside  the  sphere,  i.e.  ,  p  <1  need  be  considered. 
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Project  the  line  segment  QQ1  (see  Figure  7)  onto  the  z-axis.  The  end 

points  of  this  projection  are  Q  and  Q  '  ,  the  Z  components  respectively 
^  3  3 

of  Q  and  Q1  ,  Q  and  Q  '  form  a  trichotomy  of  the  diameter  UV.  Consider 

O  O 

the  following  cases: 


(1) 

Q  ' 

3 

< 

<  P 

(2) 

Q  1 

3 

< 

S 

=  P 

(3) 

Q  ' 

3 

< 

p  < 

Q 

3 

(4) 

Q  ' 

^3 

= 

p  < 

Q 

(5) 

p  < 

Q. 

'  < 
3 

Q 

3 

We  find,  using  Equations  (23), 


-  -  /  2 
Q3  =  Q,k=dc  +  mvl  -  d 


and 


*  r  2 

Q^'  =  Q'  •  k=dc  -  m  V 1  -  d 


(27) 


If  Q  '  ^  p  <  Q  ,  we  have 
o  o 


-  m 


J 1  -d^  +  dc— p^m  \/l  -  d^  +  dc 


then 


-  m  N 1  -  d  ^  p  -  dc 


—  m  sfl  - 


or 


(p  -  dc)  —  m 


(1  -  d2) 
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then 


(p  -  dc)2  -  m2  (l  -  d2)  ^  0  . 


(28) 


By  expanding  and  rearranging  Equation  (28)  is  equivalent  to 

F  =  (d  -  pc)2  -  m2  (l  -  p2)  ^  0  .  (29) 

Since  each  of  the  above  steps  is  reversible,  we  have  proved  the  theorem: 

2  2 

If  (p  -  dc)  ^  0  and  p  <1,  then  B  -  4AC  —  0  if  and  only  if 

2 

Q  '  —  p  —  Q  .  If  the  equality  holds,  then  B  -  4AC  =  0. 

We  have  taken  care  of  Cases  2,  3,  and  4.  We  now  consider  Cases  1 

and  5. 

Factoring  Equation  (28)  we  have 

F  =  (p  -  dc  +  m  */l  -  d2)  (p  -  dc  -  m  n/i  -  d2)  (30) 

and  from  Equations  (27) 

F  =  (p-Q3')(p-Q3).  (31) 

Now  F  >  0  if  both  factors  of  the  right  member  have  the  same  sign, 

i.  e.  , 

P  >  Q3‘  and  p  >  Q3  (32) 

or 

p  <  Q  '  and  p  <  Q  .  (33) 

o  3 
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Inequalities  (32)  give  Q  '  <  Q  <  p  ,  and  inequalities  (33)  give  p  <  Q  '  <  Q 

O  O  O  O 

and  we  have  the  theorem: 

2 

If  (p  -  dc)  *  0,  then  B  -  4AC  <  0,  if  and  only  if 
P  >  Q3  or  p  <  Q3'  . 

2 

If  B  -  4AC  <  0  and  A  *  o,  the  locus  is  an  ellipse  —  if,  in  addition, 
AA  (or  CA)  <  o.  There  is  no  locus  if  AA  (or  CA)  >  o. 

Since  A  <  o  [see  Equation  (18)],  it  remains  only  to  show  that  A  >  0, 
(or  C  >  0). 

It  has  already  been  shown  that 

B“  -  4AC  <0  if  (d  -  pc)2  +  (a2  +  b2)  (p2  -  l)  >  0  .  (34) 


Equation  (34)  leads  to 


(p2  _  i)  >  "  (d  ~  Pc)" 
VP  '  ^2  , 2'  ’ 


la“  +  b") 


then  multiplying  by  A? 


2/  2  \  -  A2i<d  -  pc>2  Ai  =  a'  b 

AjvP  V  >  /  2  2\ 

\a  +  b  /  A.  +  0  for  some  i 


and  finally, 


2  2/2  \ 

(d  -  pc/  +  A.(p  -  1)  >  (d  - 


pc) 


.  2  ,2  2\ 

2<a  +  b  -  Aj\  A.  -  a,  b 


2  2 
a  +  b 


>  0 


A.  +  0  for  some  i  . 
i 


(35) 


23 


Now, 


(d  -  pc)2  +  A2(p2  -  lj  =  A  f  for  = 


=  C  \  for  A.  =  b 
1 


and  by  Equation  (35)  A  >  0 

2 

and  C  >  0  when  B  -  4AC  <  0 
Hence  A  A  <0  and  CA  <0  and  the  locus  is  an  ellipse. 

Now  consider  the  case  |cj  =1.  If  c  =  +  1 ,  then  a  =  b  = 
Equation  (19)  becomes 

(d  -  pc)2(x'2  +  y' ")  +  [  (d  ±  l)(p  -  w)  ]  [  (d  +  l)(p  -  w)  ] 

.2/  .2  .2 


(d  -  pc)  (X1  +  y' 


(x1  “  +  y'2)  =  (l  -  d2)(p  -  w)2 


from  Equation  17 


and 


B2  -  4AC  =  -  4(d  -  pc)4 


2  4/ 

A  =  -  4(p  -  w)  (d  -  pc)  ^1  -  d  /  . 


If  (d  -  pc)  =  0 ,  there  is  no  locus. 

If  (d  -  pc)  #  0 ,  the  image  is  a  circle  with  center  (0 ,  0) 


R  = 


(l  -  d2)  (p  -  wf  

p  -  w 

J. 

(d  -  pc)2 

(d  -  pc) 

a2 . 


The  circle  is  a  null  circle  if  p  =  w  or  if  d  =1. 


0. 


=  0 


and  radius 
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SECTION  VI 


APPLICATIONS 


In  order  to  illustrate  the  above  theory,  the  projection  of  a  given 
spherical  zone  boundary  will  be  shown  for  various  values  of  p  and  w  . 

To  define  the  zone  boundary,  take 

A0  4  (i  +  j  W2  k)  and  d  =  ~ 

A  further  requirement  is  that  the  315°  meridian  lie  along  the  positive 
x  axis  and  that  the  north  pole  be  the  center  of  projection.  Thus 


C  =  k  and  B  =  (cos  315°)i  +  (sin  315°)j  =  ~  (i  -  j) 


From  Equation  (10)  we  get 


a  =  225 


and  from  Equation  (11) 


R  = 


/W2  42 

2  "  2 

42 

2  2 


0  0  1 


(37) 


Then, 


42  -  - 

A  =  *  •  A0=0  -  i+—  j  +  -k 


(38) 
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The  geometry  is  shown  in  Figure  8. 
Consider  five  cases: 


(1)  P 

(2)  P 

(3)  P 

(4)  P 

(5)  P 


0;  w  =  1, 

2»  w  =  1, 

1;  w  =  0,  and 
2;  w  =  0 


Case  1.  p  =  -  1;  w  =  0 


This  is  a  stereographic  projection.  Equation  (19)  reduces  to 

9  9 

x'  +  y'  -  y'  =0 


or 


x'2  + 


Thus  the  image  is  a  circle  with  center  at 
With 


(y'  -  I)  =  4  ' 

M) 


and  radius  — . 

2 


\[2  r  sl2  r 
A'  =  j  +  —  k 


~  _\f2  -  \[2  - 

N'  =  — ^  j  +  ~ ~  k  (see  Equation  22), 


and  finally 


-  */2  -  \[2  - 
A'  +—N’ 


=  k  , 


Q'  =  j  ; 
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A  -19,246 


(  a)  CASE  I 


(  b)  CASE  2 


(c  )  CASE  3 


(  d)  CASE  A 


IMAGE 


(e)  CASE  5 


Figure  8.  Projection  of  Spherical  Zone  Boundary  for  Various  Values  of  P  and  W 
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then  Q  =  1  and  Q  '  =  0.  Since  P  =  -  1  <  Q  '  ,  the  result  is  as  expected 

o  o  o 

from  Equation  (34)  and  Table  II. 

Case  2.  p  =  0 ,  w  =  1 

This  is  the  gnomonic  projection.  Using  a  similar  procedure  as  in 
Case  1  we  find,  since  p  =  Q  ',  that  the  image  is  the  parabola 

O 

x'2  =  2y'  • 

Case  3.  p  =  w  =  1  . 

This  is  a  "geometric"  projection.  Since  Q  '  <  p  <  Q  ,  we  find  the 

O  «J 

image  to  be  the  hyperbola 

x'2  -  2y'2  -  2y'  =  0 


2  4 


Case  4.  p  =  1,  w  =  0 

In  this  case,  p  =  Q  and,  since  the  generating  plane  contains  the 

O 

point  of  projection,  p  -  dc  =  0 .  Using  Equation  (19),  the  image  is  the 
line  y'  =1. 

Case  5.  p  =  2,  w  =  0 

Here,  p  >  Q  and  the  image  is  the  ellipse 

O 


4 
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SECTION  VII 


DISCUSSION  OF  THE  CONIC  SECTIONS 


THE  NONDEGENERATE  ELLIPSE,  HYPERBOLA,  AND  PARABOLA 

2 

For  this  case  A  +  0  and  B"  -  4AC  +  0. 

Simplify  Equation  (19)  by  using  the  rotation  given  in  Equation  (20).  This 
rotation  removes  the  x'y1  and  y'  terms.  The  resulting  equation  is 

2  T  2  2/  2  M  2  2 

x'  |_(d  -  pc)  +  m  (p"  -  ljj  +  y'  (d  -  pc)"  -  2mx'  (p  -  w)(pd  -  c) 


(  2  2\  2 
+  (d  -  c  /  (p  -  w)  =  0  . 


(39) 


Noting  that 

(d  -  pc)"  +  m“(p2  -  l)  =  (pd  -  c)2  +  (l  -  d2)(p2  -  l) 
put  Equation  (39)  in  a  standard  form, 


m(p  -  w)  (pd  -  c 

■1 

(pd  -  c)2  +  1 

(i-a2) 

to 

(p  -  w)"(l  -  d")  (d  -  pc)~ 


(P  ~  w) 


2  (l  -  d2) 


=  1 


[(Pd  -  c)2  +  (l  -  d2)(p2  -  l)]"  (pd  -  c)2  +  (l  -  d“)fp2  -  j) 


(40) 


Now,  with 


F  =  (pd  -  c)2  +  (l  -  d2)(p2  -  l)  =  (d  -  pc)"  m" (p2 
=  (p  -  dcf  +  m2(d2  -  l)  , 


\ 

V 
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the  center  of  the  conic  is 


H<V 


y  ' )  where, 
0 


,  _  m(p  -  w)  (pd  -  c) 

A  A  __ 

0  F 


and 


The  axes  of  the  conic  sections  given  by  Equation  (40)  are  such  that  one 
axis  is  always  on  a  radial  line  from  the  center  of  projection.  This  is 
referred  to  as  the  "radial  axis, "  and  the  length  of  the  semiradial  axis  is 


r 


=  4l-z  y}. 


(l  -  d2) 


(41) 


The  axis  normal  to  r  is  the  "transverse"  axis,  and  the  length  of  the  semi- 
transverse  axis  is 


Note  here,  that  when  F  <  0,  the  conic  is  a  hyperbola,  and  for  F  >  0, 
the  conic  is  an  ellipse. 


Using  the  inverse  rotation  to  Equation  (20), rotate  the  system  to  its 
original  position.  The  position  vector  of  the  center  of  the  conic  thus  is 


=  a(P  -  w)(pd  -  c)  -  b(p  -  w),(pd  -  c)  -  c 

P  P  J 


(43) 
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and,  in  the  image  plane,  the  cartesian  coordinates  of  the  center  are 


x 


,  _  a(p  -  w)  (pd  -  c) 


H 


and 


,  _  b(p  -  w)  (pd  -  c) 


H 


(44) 


Since  c  is  invariant  under  the  rotation  used,  the  semiradial  and 
semitransverse  axes  are  as  given  by  Equations  (41)  and  (42),  respectively. 

F  <  0 

If  F  <  0,  the  conic  is  a  hyperbola. 

The  coordinates  of  the  foci  f  and  f  are 

X  t-d 


and 


f  =  x  '  +  S 
l,x  H  m 


U  „  =  y 


f  =  X  '  -  Si 

2 ,  x  H  m 


Sk 


1 ,  y  H  m 


(45) 


f0  =  y  '  -  —  for  m  *  0  ; 
2 ,  y  H  m 


(46) 


and 


g 


_  /( P  ~  wf(l  -  d")(d  -  pc)  (P  -  w) 


2  ■  -2(i  -  a2) 


ELr-»  (i  _  d2)  Jl  I  p2  ; 

r 


note:  -  1  <  p  <  +  1 


(47) 
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F  >  0 


If  F  >  0,  the  conic  is  an  ellipse. 

If  r  >  t,  the  coordinates  of  the  foci  f  and  f  are  the  same  as  given 

i.  m 

by  Equations  (45)  and  (46)  except  that 


S  = 


(p  -  w y 


2(l  -  d2) 


2  ■  -2(l-d2) 


(d  -  pc)  -  (p  -  w) 


(48) 


If  r  <  t,  the  foci  are  on  the  transverse  axis.  To  locate  these  foci 
consider  the  ellipse  as  given  in  Equation  (40). 

The  position  vector  of  the  center,  H,  is 

H'  =  Xq1  i  +  wk  . 

If  x  '  >  0  ,  the  focus  which  is  angularly  displaced  from  the  radial  axis 
in  a  positive  direction  has  the  position  vector 

fl  +  "  x0'  1  +  gJ  +  ’ 

and  the  other  focus  is  given  by 

f2  _  =  V  i-gj  +  wk- 

If  x  '  <  0,  the  direction  of  the  angular  displacements  are  reversed. 

Using  the  inverse  rotation  as  in  Equation  (20),  rotate  the  system  to  its 
original  position  and  find 

V  +  =  m  [(“o'  '  bg)  i+(bxo’  +  "*)  ?! 
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and 


(50) 


f 

2, 


-  =  m  [(“o'  +  bg)‘+(bxo' 


Since  r  >  0,  and  t  >  0,  to  determine  the  conditions  for  which  r  <  t,  let 


or 


then 


< 


(d  -  pc) 


2 


<  F  =  (d  -  pc) 


2 


+ 


2/2 

m  (p 


and  we  have 

2/  2  \ 

m  (p  -  lj  >  0  .  (51) 

2 

Thus  for  m  =£  0,  p  >  1  is  the  condition  for  which  the  transverse  axis  is  the 
longer  axis. 

If  r  =  t,  the  conic  is  a  circle.  The  conditions  for  which  r  =  t 
follow  from  Equation  (51); 

m  (p  -  l)  =  (1  -  c  )(p  -  l)  =  0  , 


gives 


c  =  ±  1 
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and 


P  =  ±  1 


(52) 


as  the  required  conditions. 

Thus  a  circle  is  possible  only  in  the  stereographic  projection  (p  =  ±  1) 

A  A 

or  if  A  =  k,  i.  e.  ,  if  the  generating  plane  is  perpendicular  to  the  z-axis. 

The  radius  of  the  circle  is 

,,  _  |(p  -w) (d  -  pc)  |v'l  -  d2 

|F| 

For  F  =  0  and  A  ^  o,  the  conic  is  a  parabola. 

Using  the  rotated  system,  Equation  (39)  reduces  to 

2  2  /  2  2\  2 

(d  -  pc)  y'  -  2mx'  (p  -  w)(pd  -  c)  +  Vd  -  c  J(p  -  w)  =0  (53) 


or  in  a  standard  form, 


1 2  _  2m(p  -  w)  (pd  -  c)  , 

y  -  2 

(d  -  pc)  L 


2V  x2 
-  c  j(p  -  w) 

2m  (pd  -  c) 


This  parabola  has  a  center,  H(x^' 


y  '  ),  with 
J0 


x 


0 


,d"  -  c2)  (p  -  w) 

2m  (pd  -  c) 


(54) 


(55) 


and 
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The  coordinates  of  the  focus  are 


(x0'  +  h,  0)  (56) 

where 

h  =  m(p  -  w»pd  -  c>  ;  note:  JI  =  focal  length.  (57) 

2(d  -  pc) 

The  center  of  the  parabola  in  this  rotated  system  is  given  by 

H'  =  x  '  i  +  0  •  j  +  wk  . 

0  J 

Rotating  the  system  to  its  original  position,  as  before,  we  find 

H  =  R"1  H  1 1  =  x  '  ~  i  +  x  '  —  j  +  wk  .  (58) 

0  m  0  m  '  1 


The  cartesian  coordinates  of  the  center  (in  the  image  plane)  are 


x 


H 


(2  2\ 
ayd  -  c  )  (p  -  w) 

2 

2m  (pd  -  c) 


and 


2 

2m  (pd  -  c) 


(59) 


Using  a  similar  procedure,  the  coordinates  of  the  focus  are 


x 


f 


m 


a 
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and 


(60) 


THE  DEGENERATE  CASES,  A  =  0 


The  following  cases  are  considered. 


Case  1.  p  =  w. 


In  this  study  p  #  w. 
Case  2.  (l  -  d^)  =  0. 


Project  a  point  ellipse  and  if  c  +  1  the  image  is  given  by  Equation  (44). 

If  c  =  +  1,  we  have  no  image  if  (d  -  pc)  =  0;  otherwise,  the  image 
is  the  point  (0,  0). 

Case  3.  d  -  pc  =  0 

Using  the  rotated  system,  Equation  (39)  reduces  to 
x'2  m2(p2  -  l)  -  2mx'  (p  -  w)(pd  -  c)  +  (d2  -  c2)  (p  -  w)“  =  0  (61) 

Solving  Equation  (61)  we  have 


x 


(P  -  w)(pd  -  c) 

m(p2  -1) 


for  p  #  1  . 


(62) 


2 

If,  in  addition  p  =1,  Equation  (61)  reduces  to  (for  m  +  0) 


x 


(63) 
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and 


x 


,  (d  -  c)(l  +  w) 

=  - ^ -  for  p  =  -  1 


(64) 


To  find  the  equation  of  the  image  in  the  original  system,  let  d^  =  di  +  wk 
be  a  vector  normal  to  the  lines  given  in  Equation  (62)  and  Equation  (63),  where 
d  #  0  is  the  directed  distance  from  (0,  0,  w)  to  the  line. 

Then,  rotating  d  into  the  original  system  we  have 


-r  Da:  Db  ?  r 
d  =  —  i  +  —  i  +  wk. 
m  m 


The  projection  of  d  on  the  x'  y'  plane  is 


TT  Da  r  Db  t 
d  =  —  l  +  —  i 
m  m 


then  x  •  ^  =  D  gives  the  equation  of  the  line  as 


+  —  y'  =  D  . 
m  m 


The  images  corresponding  to  Equations  (62),  (63),  and  (64)  are 
for  Equation  (62) 


,  .  ,  ,  (P  -  w)(pd  -  c) 

ax'  +  by'  =  ,  x - L~  ’ 

(P  -  1) 


for  Equation  (63) 


ax'  +  by'  -  g-LPH1  -  w>  for  p  =  +  1. 
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and  for  Equation  (64) 

ax'  +  by'  =  ^-r-£KLlw)  for  p  =  .  ! 

Case  4.  c  =  d  =  0. 

In  this  case  the  zone  boundary  and  the  point  of  projection  lie  in  the 
plane  ax  +  by  =  0 .  The  image  obtained  is,  then,  the  intersection  of  this 
plane  with  the  image  plane,  Z  =  w.  The  image  is,  then,  the  line 

ax'  +  by'  =  0  . 
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SECTION  VIII 


POSITIONING  THE  IMAGE  PLANE 


There  is  no  restriction  on  the  analysis  if  p  —  -  1. 

A  point  T(x,  y,  z)  of  the  sphere  with  position  vector  T  =  xi  +  yj  +  zk 
has  a  projection  on  the  (x,  y)  plane  given  by 


T 

xy 


=  xi  +  yj 


The  position  of  the  image  of  point  T  is 


T' 


p  -  w 
p  -  z 


(xi  +  yj)  . 


Thus,  if  (p  -  w/p  -  z)  >  0,  the  image  will  have  the  same  direction 
from  the  origin  as  the  preimage,  otherwise  the  direction  will  be  increased  by 
180  degrees. 

The  factor  (p  -  w/p  -  z)  >  0  if  and  only  if 


w  >  p  when  z  >  p 


or 

w  <  p  when  z  <  p  . 

Now,  if  -  1  —  p  <  Q  '  then  all  points  (x,  y,  z)  of  the  preimage  have 

O 

z  >  p  hence  we  require  w  >  p . 

If  Q  1  <  p  <  Q  ,  then,  for  some  points  of  the  preimage,  z  <  p,  and 
3  3 

only  that  portion  of  the  preimage  for  which  z  >  p  shall  be  used.  We  require, 
al  so ,  that  w  >  p  . 
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If  p  >  1  >  Q  ,  then  for  all  points  (x,  y,  z)  of  the  preimage,  z  <  p, 

O 

and  we  require  then  that  w  <  p  . 

If  p  >  1,  project  only  that  portion  of  the  preimage  for  which  z  ^  1/p . 
It  can  be  shown  that  a  right  circular  cone  whose  apex  is  at  P  is  tangent  to 
the  sphere  on  the  circle  given  by 


and 

1 

z  =  —  . 
P 


We  shall  project  those  points  of  the  preimage  for  which  z  >  l/p  and 
we  shall  require  w  <  p  . 
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APPENDIX 


DERIVATION  OF  THE  DISCRIMINANT  A 

The  expansion  of  the  determinant  giving  A  is  a  rather  tricky  algebraic 
exercise.  For  this  reason  the  essentials  are  given  below. 

Substituting  from  Equation  15 

2  jjd- pc)2+ a2(p-  l)j  2ab(p2-l)  2a(c  -  pd)  (p -  w) 

2ab(p2  -  1 )  2  [jpc  -  d)“  +  b2(  p“  -  l)]  2b(c  -  pd)  (p  -  w) 


A  =  1/2 


2a  (c  -  pd)  (p  -  w) 


2b(c  -  pd)  (p  -  w) 


(d2-c2) 


(p  -  W) 


An  obvious  factoring  gives, 


A  =  4  (p  -  w) 


(d  -  pc)2  +  a^(p~  -  l)  ab(p2  -  l)  a(c  -  pd) 

ab(p“  -  l)  (pc  -  d)“  +b2(p2  -  l)  b(c  -  pd) 

a(c  -  pd)  b(c  -  pd)  (d2  -  c2) 

2  2 

Factoring  out  a  b  ,  then  subtracting  column  2  from  column  1.  then  factoring 

2 

out  (d  -  pc)  ,  and  subtracting  row  1  from  row  2,  we  have 


2  2  2  2 
A  =  4a  b  (p  -  w)  (d  -  pc) 


A  .  (p2-/ 


J_  +  M  M  -  pc) 

2  .2  i  .2 

a  b  /  b 


(c  -  pd) 


c  -  pd 
0 


,2  2 

d  -  c 
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Finally,  expanding  on  the  first  column, 


A  =  -  4(d  -  pc)  (p  -  w)‘ 


(l  -  o2)- 
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